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Abstract 

Spatially homogeneous solutions of the Landau-Lifshitz-Gilbert equation 
are analysed. The conservative as well as the dissipative case is considered ex- 
plicitly. For the linearly polarized driven Hamiltonian system we apply canon- 
ical perturbation theory to uncover the main resonances as well as the global 
phase space structure. In the case of circularly polarized driven dissipative 
motion we present the complete bifurcation diagram including bifurcations up 
to codimension three. 
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1 Introduction 



The investigation of strongly driven ferromagnetic systems is an interesting field of 
research especially from the point of view of nonlinear dynamics. The fundamental 
equation of motion, which governs the dynamics of a macroscopic magnetization has 
been proposed already by Landau and Lifshitz [0] 

S = -S X (He//(S, t) + rS X He//(S, t)) . (1) 

It describes the motion of a macroscopic magnetization S under the influence of 
an effective local magnetic field He//(S,t). The second term in this equation, fre- 
quently called the Gilbert damping, represents the dominant contribution among 
several dissipative terms, which arise in microscopic derivations of this equation 
(e.g. ref.@, It is a special feature of eq.(|l|) that it preserves the modulus of 
the magnetization. In order to give eq.(||) a definite meaning one has to specify 
the effective field. Usually contributions from external fields, dipolar and exchange 
interactions, anisotropy etc. are taken into account. They turn eq.(0) in general 
into a complicated partial differential equation which may even be nonlocal in space. 
Hence it is impossible to discuss the dynamical behaviour in general. We will focus 
in our article on the analysis of spatially homogeneous states. Even this seemingly 
simple problem turns out to be highly nontrivial. 

To begin with the effective field has to be specified. We take a static external 
field, a time dependent transversal driving field and an uniaxial anisotropy parallel 
to the static field into account while contributions from the exchange interaction 
obviously vanish and the dipolar interaction may be thought to be contained in 
the anisotropy term |Q. Under these conditions eq.(|l]) becomes a two dimensional 
system of nonautonomous ordinary differential equations, which is Hamiltonian or 
dissipative depending on the value of the damping constant. We will treat both 
cases in our paper separately. 

Discussions of low dimensional Hamiltonian spin systems can be found in the 
literature. A treatment of undriven systems may be found in |^. Furthermore 
Frahm and Mikeska have shown that the equations of motion are integrable, if 
a circularly polarized driving field is applied to a single spin. For linearly polarized 
driving fields their numerical simulations indicate that the system is not integrable. 
In contrast to their approach we will analyse the Hamiltonian dynamics using canon- 
ical perturbation theory and compare the results with numerical simulations. 

On the other hand a systematic analysis of driven dissipative spin systems does 
not seem to be available in the literature. We will study the damped spin system 
under the influence of a rotating driving field. In that case the system becomes au- 
tonomous using a transformation into a rotating coordinate system. Although such 
a system cannot become chaotic we will show by performing a careful bifurcation 
analysis that highly nontrivial phenomena occur. 
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2 Hamiltonian Dynamics under Linearly Polar- 
ized Pumping 



Undamped spin systems provide a nontrivial example of simple Hamiltonian systems 
0. It is well known that an axisymmetric system is integrable if a rotating driving 
field is applied. Hence we focus in this section on systems which are driven by linearly 
polarized fields. They give rise to very complicated dynamics. To be definite the 
effective field is chosen as 

Heff{S,t) = [Hz - aSz)ez + h±cos{ujt)e^ . (2) 

Spin systems of this type have been discussed already in the literature. But in 
contrast to former approaches 0, § a static field is incorporated in our treatment 
which breaks the plane symmetry Sz — > —Sz- For that reason the methods employed 
in these references cannot be applied to our equations of motion. 

We will use canonical perturbation theory for weak driving fields to uncover 
the phase space structure of our system. The calculations are then confirmed by 
numerical simulations. 



2.1 Canonical Perturbation Theory 

The dynamical system (|l]) subjected to the effective field @) is described by the 
Hamiltonian 

= 7^0 + £ "Hi 

= Sz [Hz - ^) + £ /i± - SI cos(^) sin(cut) . (3) 

Here canonical coordinates S^+iSy = ^1 — S'^ exp(iy9) and the formal expansion pa- 
rameter e have been introduced for convenience. The unperturbed i.e. the undriven 
system, e = 0, is obviously integrable and admits the solution 

Si'\t) = Si'\0) (4) 
(^W(t) = y.(°)(0) + ^](^f (0))t (5) 
n{Sz) := {H,-aSz) (6) 

corresponding to the precession of a single spin. 

It is the scope of canonical perturbation theory to construct a sequence of trans- 
formations which map the original Hamiltonian approximately to an integrable one. 
Formally this is achieved by eliminating fast variables (in our case ip and t) up to a 
certain order e"^ in the Hamiltonain . In each order of this approximation finitely 
many small denominators will appear, so that on a finite collection of surfaces the 
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new Hamiltonian is not well defined. Outside a small neighbourhood of these sur- 
faces the transformed system is integrable and may serve as a good approximation 
to the original one. The small denominators signal the resonances in the system and 
have to be treated separately. They prevent the expansion from being convergent 
and even from being valid asymptotically. 

We perform the perturbation theory by using a modified version of the von 
Zeipel method. The perturbation expansion is formalized with the help of the Lie 
technique^]. We are looking for a canonical transformation 5*2, Sz,if which 

transforms the Hamiltonian @j into an integrable one. Obviously the transformation 
depends on the expansion parameter e. If w{J,ip,e) denotes the generator of the 
infinitesimal transformation it obeys 

f^ = [xM . (7) 

where x stands for and (p, respectively and [., .] denotes the Poisson bracket. 
The evolution operator x{e) = Tx of this formal Hamiltonian system with "time" 
e defines the desired transformation. With the help of the Lie operator L = [w, ] 
eq.(0) can be cast into the form 

dT 

This algebraic version of the transformation is well suited for perturbation expan- 
sions. Using the formal power series 

oo 

L = $]£"L„+i (9) 

n=0 

oo 

T = J^e^T^ . (10) 

n=0 



eq.(H) can be solved for T 



n — 1 

n 



m=0 



— X! Ln-mT„^ . (12) 

m=0 



1-1 T rf-i-l 

n 



Of course the generator w has to be specified in order to give these expressions 
an explicit meaning. This is accomplished by considering at first the transformed 
Hamiltonian 

n{x{x,e)) =H{x) (13) 

good introduction to the field of Hamiltonian perturbation theory is given in [|| chapter 
2 and [1^ chapter 5. For a detailed introduction to the Lie transformation methods see e.g. |^ 
section 2.5. 
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which in terms of the evolution operator (|l^) is given by H = T ^7i. From the 
formal series expansions 

oo 

H = Y.^'^'^n (14) 
n=0 

oo 

H = Y.^^'^n (15) 
n=0 

one obtains the recurrence relations 

dw, 



n-l 



dt 



m=l 



+ [Wn, Uo] = n [Hn - Hn) - Y i^n-m'Hm + mT^^^Hm) ■ (16) 



Owing to the fact, that the unperturbed motion can be integrated explicitly (cf. 
eqs.(i[^) these equations can be solved successively once Hm is fixed in an appropri- 
ate way. It is this freedom which is used in the sequel to construct the perturbation 
series. 



2.1.1 The First Order 

Let us start with the simple first order perturbation expansion. Eq.(^) reads ex- 
plicitly for n = 1 

^ + [wuno] = Hi -Hi 



= ni-h±^l- SI cos(^) sin(cut) . (17) 

In order to render the new system integrable one requires that the new Hamiltonian 
depends solely on the action variable. Observing this condition Tii is chosen in 



such a way that the solution of eq. (|l7|) possesses no secular contributions. This is 



achieved by choosing Tii = 0. We find for the generator of the transformation 

' cos {i-p + ujt) COS ((^ — cut) 



wi = h±Jl - SI — / (18) 

and hence for the transformation itself 

Tf=l-e[wiJ] . (19) 

It has already been mentioned that the perturbation expansion leads to small de- 
nominators indicating resonances. They appear at ^{Sz) = ±uj. Higher harmonics 
do not occur in the present case because of the special nature of the coupling of the 
driving field. 
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2.1.2 Removal of First Order Resonances 



The local structure of the phase space near the resonances cannot be described by 
the approach of the previous paragraph. However we can uncover the local phase 
space structure near some special resonance by using a different choice of the first 
order Hamiltonian T^i. To be definite let us consider the resonance at Q = —uo. In 
order to motivate the new choice we transform into a rotating frame0 using the new 
variables = Sz and (p = ip + ut. Then the Hamiltonian (|^) reads 



n = Sz \^Hz - ^ + ujj + e h^^Jl - S^^ (sin (^) + sin {ip - 2 ut)) . (20) 

Now the new angle <^ is a slow variable near the resonance under consideration. 
Thus we do not try to remove it from the Hamiltonian. We are now looking for a 
solution of eq.(|l6l) which in the present case reads 

dwi 
Choosing 



+ [wu Ho] = ^1 - /i± V 1 - S^z (sin (0) + sin {0-2 ut)) . (21) 



ni = h^^l-S^, sinip (22) 

one ensures that the solution 

cos {ip — 2 ut) 



2 IHz — aSz — uj I 



is regular near the resonance. In addition the new Hamiltonian reads 



n = Sz\Hz-^ + ujj+e h^^Jl - SlsiYi0 + O [e') (24) 

is independent of time up to the first order in the expansion parameter and represents 

an integrable system. The phase space structure of this system is sketched in Fig.l. < Fig.l 

The system has a hyperbolic fixed point dX (p = and an elliptic fixed point at 

(p = ^. In order to describe the motion near the resonance 5*^ = {Hz — uj)/a, we 

expand the Hamiltonian in powers of small deviations ASz from the resonance. In 

lowest order one obtains the well known Hamiltonian of a pendulum 

H^ + u a/,A\2 , I /Hz - cj\^ , , 
2^ V \~^^) "'""^ ■ ^^^^ 

This expression describes the generic case of near resonance motion and is studied 
in the literature (e.g. ref. [Il|) very well. 



^It should be mentioned that the introduction of the rotating frame is by no means necessary. 
However the apphcation of these coordinates is physically more appealing. 
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2.1.3 Second Order Islands 



Taking the higher order contributions of the Hamiltonian (^) into account leads 
to perturbation terms of higher order in e and ASz- Especially a time dependence 
of frequency 2uj will enter the Hamiltonian. We note that the internal frequency of 
system ( ^5|) is of the order O (i/e), whereas the external frequency is of order O (1). 
Hence the full system will develop resonances of order where the integer p is of 
the order of magnitude e~^^'^ . Now these second order resonances themselves result 
in p hyperbolic fixed points and the heteroclinic tangles belonging to them. The 
amplitudes of these resonances is of the order O {l/ (^e"^/^^!^. For a much more 
detailed analysis of this case we refer the reader to the literature (e.g. [§] section 
2.4b). 



2.1.4 The Second Order 

The previous paragraphs have dealt with the analysis of the first order resonances 
extensively. If we proceed to higher order additional resonances will occur. Following 
the lines of paragraph 2.1.1 we obtain in second order from eq.flTB) 

+ [w2,no] = - [wi,ni] 

(cos i2cp) + cos (2^t) - 1) (1 - S',) jn' + uj') 

= A /To — il I O n 

^ (2 (1]2 - cj2))2 



+h 



2 Sz{l-cos{2iut))n 



The choice 



(2 (n^-iu^)) 

(l _ s^) ( ^osi2^ + 2.t) cos(2^-2^t) \ _ 

^ 8(fi2-u;2)2 ^4(fi2-cu2) ^ ^ 



ensures that the solution possesses no secular terms. The generator becomes 

,^nS,sm(2ujt) ,2 (1- S^^(n^ + uj^)sm(2ip) 
W2 = -h\- — -— — + h\a- 



+hi a- 



2 {I- SDiSl'^ + oj^)sin{2ujt) 



8uj (fi2 _ 

Clearly an additional resonant small denominator appears at Vt[Sz) = 0. The full 
transformation reads up to the second order in view of eq. (p!T|) 

Tf = l-e [wu f]-j i[w2, f] - [^1, [^1, /]]) . (29) 
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and reflects the first and second order resonances. 

In analogy to the procedure described in paragraph 2.1.2 the new resonant de- 
nominator can be removed by a different choice for the new Hamiltonian. The 
appropriate expression is given by the time average of the right hand side of eq.(^) 



(cos {2ip) - 1) (1 - S^,) + uj^) 



The corresponding solution for the generator remains regular near the second order 
resonance 

+hla (l - S!) f ~ + ^""^ ^ 't^ ] ■ (31) 

It should be mentioned that one has to pay for the removal of the resonance at 
= by the explicit occurrence of the angle variable ip in the expression (|30| ) for 

Considering now the transformed Hamiltonian and recalling that Tii vanishes 
(cf. paragraph 2.1.1) we have in view of eq.(^) 

n = no + e^n2 . (32) 

Expanding this Hamiltonian near the second order resonance ^l{Sz) = in powers 
of small deviations AS^ from the resonant surface we obtain again the Hamiltonian 
of the pendulum in the lowest non vanishing order 



2 8ato'^ 



n = — {AS,y-e'hi ' „ cos2y^ . (33) 



In contrast to the previous case of the first order resonance the width of the resonant 
region scales with the amplitude of the driving field eh±. 



2.1.5 Global Analysis 

The discussion of the previous paragraphs has shown that one can detect resonances 
in the spin system and can describe the phase space structure near the resonances 
approximately. In addition it would be of course desirable to gain a global overview 
over the phase space structure. This goal however cannot be achieved with the 
help of the local analysis of the previous sections because the generators contain 
singularities away from the resonances (cf. eq. (|23|) ). 
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In order to get an overview of the phase portrait of the Hamiltonian system 
we follow a method which has been proposed by Dunnet at al. ||l2|. Its main 



idea is quite simple and is based on the construction of an approximate constant of 
motion. Suppose that the original Hamiltonian has been transformed in some order 
to a system whose Hamiltonian depends only on the action variable Sz- Then an 
arbitrary function of the action f{Sz) yields a constant of motion of the transformed 
system. Hence the counter image T~^/ is also a constant of motion and contains 
the topology of the phase space curves in the original variables. This naive view 
is in general meaningless because of the small denominators involved. But the 
construction can be performed with the singular expression in every order of the 
perturbation theory if the constant of motion / is suitably chosen. 

To be definite consider the perturbation expansion up to second order (cf . eq. 
Formal application of expansion 

(IH) leads to 



+ 



Wi, 



(34) 



This expression contains singularities at the resonances ^l{Sz) = 0, ±uj. If however 
the function / is chosen in such a way that the zeros of its derivative cancel all these 
singularities the expression remains regular in the whole phase space. The authors 



of ref.|]12j suggested this expression to be a reasonable approximation for a global 
invariant of the system. 

In principle one can use any function / which satisfies this constraint. For 
graphical purposes it is however convenient to use an expression whose numerical 
values are of the same order of magnitude in a region of phase space as large as 
possible. We therefore start from the following expression 



cos 



2uj 



Q{x) dx 



(35) 



Combining eqs.(|l8[), ( |28[) , and (^) the invariant up to second order perturbation 
theory is obtained after some tedious calculation. We refrain from writing down the 
lengthy result explicitly. 

The level lines of the invariant constructed in this manner provide a good ap- 
proximation of the phase space structure. We give plots of these level lines at fixed 
time t which is equivalent to a Poincare section at the same time. Fig.2 shows the < Fig.2 
lines for the parameter values eh± = 0.01, a = 1, Hz = 0.5 and u = 0.2. Reso- 
nances of first and second order are clearly visible. The size of these resonances is 
in accordance with the discussion given in the preceding paragraphs. Additionally 
second order islands can be recognized near the first order resonance. But these 
structures are reproduced only qualitatively because they depend on the expansion 
parameter in a highly nonanalytical way (cf. paragraph 2.1.3 and ref. [|I^). If the 
amplitude of the driving field is increased the level lines which cross the phase space 
are destroyed. Fig. 3 shows an example for the parameter values h± = 0.025, a = 1, < Fig.3 
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Hz = 0.5 and on = 0.2. This effect may be attributed to an overlap of low order reso- 
nances and signals the destruction of KAM tori. However higher order perturbation 
calculations would be necessary to obtain more reliable results. 



2.2 Numerical Results 

We confirm the analytical perturbational calculations of the previous section by di- 
rect numerical simulations of the full Hamiltonian dynamics. For that purpose the 
system has been integrated by using the subroutine D02BAF of the NAG library. 
Poincare plots have been computed by integrating over 100 periods of the driving 
field starting from 400 initial conditions distributed uniformly at if = 0,±7r/2,7r. 
Fig.4 shows the Poincare plot for the parameter values used in Fig.2. Good agree- < Fig.4 
ment is observed. But e.g. third order resonances are found in the numerical 
simulation being not incorporated in the perturbative approach which is of second 
order only. Fig. 5 contains data at the parameter values used in Fig. 3. Obviously < Fig.5 
KAM tori are destroyed by resonance overlap. The agreement between the pertur- 
bation expansion (cf. Fig. 3) and the numerical simulation is merely qualitative. It 
signals that the validity of the expansion breaks down if the system is chaotic in 
large regions of the phase space. 



3 The Circularly Polarized Driven Dissipative Sys- 
tem 

We will now study the dynamics of the dissipative system which is fundamentally 
different from the case analysed in the previous section. It is our intention to 
give a rather complete survey over the dynamical behaviour in a simple but highly 
nontrivial situation. We consider a system being driven by a rotating transversal 
field. The effective field reads 

Heyj(S, t) = {Hz — aSz) + (cos {ut) e^, + sin {ut) By) . (36) 

The explicit time dependence of the equations of motion (0) may be eliminated 
using a transformation to a rotating frame because of the chiral symmetry of the 
dynamics. We obtain the autonomous equations of motion 

^ = -S X (U^fj + r S X (He// + ^e,)) (37) 

where 

Heff = {Hz-aSz-uj)ez + h±e^ (38) 

denotes the time independent effective field in the rotating frame. It is the physical 
nature of the driving field which is responsible for the different magnetic fields in 
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the reversible and the Gilbert term. This structure prevents eq. (|37D from possessing 
a Lyapunov function and ultimately causes the complex behaviour of its solutions. 
On the other hand the modulus of the magnetization is preserved in the course of 
time so that the phase space is the two dimensional sphere. This fact enables us to 
keep the discussion to a great extent analytical but at the same time prevents the 
system from becoming chaotic. 

As a constant factor can be absorbed into the time scale we restrict the subse- 
quent treatment to the case |a| = 1 without loss of generahty. We will start our 
analysis with a brief discussion of the fixed points. Subsequently local and global bi- 
furcations will be analysed which will result in a complete overview of the dynamics 
of the system. 



3.1 The Fixed Points 

Even the calculation of the fixed points gives rise to a system of algebraic equations 
in general which cannot be solved without resorting to numerical methods. But our 
fixed point problem (cf. eq.(|37|)) can be reduced to a single algebraic equation of 
fourth order 

(l - S!) {{aS^ - + ^'hlS!) - hlSl = . (39) 



Here the new parameters 

7 := (41) 

have been introduced and will be used in the sequel. The remaining components 
are determined via 

Sy = -l{l-S!) . (43) 

It is worthwhile to mention that because of eq.(|3^) the system has at least two and 
at most four fixed points. The explicit discussion is postponed to the next sections. 



3.2 Method of Analysis 

Before we enter the discussion of the various bifurcation diagrams it is useful to give 
a brief survey over the symmetries of the system under consideration. It is easy to 
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check that the equations of motion are invariant with respect to the following two 
transformations 



Ti : A ^ -A 7-^-7 a-^ -a -h_i_ t -t (44) 

T2: A ^ -A 7^-7 S,^-S, Sy = -Sy . (45) 

Taking these symmetries into account we can restrict the following bifurcation anal- 
ysis to the case a = — 1 and 7 > 0. In addition we will not analyse the full four 
dimensional parameter space (F, A, 7, h^). We restrict ourselves to the fixed value 
of the damping constant F = 0.1. 

Let us first describe the procedure which is applied to treat the local codimension 
one bifurcations (Hopf and saddle node bifurcations). It is convenient to perform this 
analysis by referring explicitly to the fact that the phase space is two dimensional. 



We rewrite eg. (1371) in terms of the complex variable Z = {S^ + iSy) / (1 + 5*2). This 
corresponds to the well known stereographic projection of the sphere to the plane 
(cf. ref.Q) 

f = fiZ)■.= {^-T)(^^h^Z + AZ-az\^^-^{l-\Z\^)j . (46) 

Saddle node bifurcations of the fixed points are determined by a single vanishing 
eigenvalue of the linearized system i.e. by the equations^ 

/(Z)=0, det(D/(Z))=0 . (47) 

The solutions of this set of algebraic equations define codimension one manifolds 
in parameter space on which saddle node bifurcations take place. With the help 
of the program PITCON [|l^ these manifolds can be computed easily by a contin- 
uation technique. In the same way Hopf bifurcations are detected. Owing to the 
dimensionality of the phase space they are determined by the equations 

f{Z) = 0, Tr(D/(Z)) = . (48) 

In addition using a well known formula for the coefficient of the third order in the 
Hopf normal form [|l^ we decide whether the bifurcation is sub- or supercritical. 
We stress that our explicit knowledge of the number of fixed points enables us to 
detect all local codimension one bifurcations. 

The discussion of the local codimension one bifurcations will show that our sys- 
tem has at most one saddle point. Homoclinic bifurcations occurring eventually in 
connection with this point have been investigated using a numerical computation of 
the distance between its stable and unstable manifolds. This method seems to be 
more suitable than the quest for limit cycles as suggested by Doedel (ref. [|T3)- The 
respective codimension one bifurcation sets are again computed via PITCON. 

^Inspection of the polynomial eq.(|39|) shows that this condition coincides with a degeneracy of 
its zeroes. 
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Finally we have concentrated on the investigation of saddle node bifurcations 
of limit cycles. It is well known that the proper bifurcation manifold is born in a 
degenerated Hopf bifurcation. Using the program AUTO we have been able to 
compute the saddle node bifurcation manifold of limit cycles also. 



3.3 Bifurcation Diagrams 



The results of our investigations are summarized in Figs. 6, 13-16. In order to get 

insight into the structure of the full three dimensional parameter space we show 

cross sections for several values of /i^. Let us first describe in detail the situation 

for hx^ = 0.1 (cf. Fig.6) . Apart from the codimension one manifolds we find the < Fig. 6 

following bifurcations^ of codimension two: 



• two cusp points at Ci and C2 

• one degenerated Hopf bifurcation at H where a saddle node bifurcation for 
limit cycles meets a Hopf bifurcation 

• two Arnold-Takens-Bogdanov bifurcations at Ai and A2 where saddle node, 
Hopf, and homoclinic bifurcations meet 

• two saddle node connections at 5*1 and S2 

• one degenerated homoclinic bifurcation^ at B where a saddle node bifurcation 
line of limit cycles meets a homoclinic bifurcation. The lines to the left of B 
cannot be resolved on the scale of Fig. 6 (cf. the small insert). 



The local and global codimension one lines divide the parameter space in ten regions. 
The typical dynamical behaviour in each region and on the boundaries is described 
in the sequel: 



• In region 1 there exist only one stable and one unstable fixed point. 

- The saddle node bifurcation line, which separates region I from region II, 
generates a saddle point and a stable fixed point. 

• Region II contains four fixed points, one saddle point, one unstable fixed point 
and two stable fixed points. 

- The saddle node bifurcation line, which separates region II from region HI, 
destroys a stable fixed point and the saddle point. 

• In region HI there exist one stable and one unstable fixed point. 

^For a detailed explanation of these basic bifurcations the reader is referred to [0 
^For a description of this bifurcation see e.g. [p^ . 
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The subcritical Hopf bifurcation line, which separates region III from region 

IV, generates an unstable limit cycle. The unstable fixed point becomes stable. 

region IV contains two stable fixed points and one unstable limit cycle. 

The subcritical Hopf bifurcation line, which separates region IV from region I 
destroys the unstable limit cycle. One stable fixed point becomes unstable. 

The supercritical Hopf bifurcation line, which separates region IV from region 

V, generates a stable limit cycle and turns one stable fixed point into an 
unstable fixed point. 

In region V there exist a stable and an unstable fixed point and a stable and 
an unstable limit cycle. 

The hmit cycles just mentioned are destroyed at the saddle node bifurcation 
line for limit cycles, which separates region V from region I. 

The saddle node bifurcation line, which separates region V from region VI (cf. 
the small insert in Fig.6), generates a stable fixed point and a saddle point. 
Moreover this bifurcation destroys the stable limit cycle. 

Region VI contains four fixed points (one saddle, one unstable and two stable 
fixed points) and an unstable limit cycle. 

The homoclinic bifurcation line, which separates region VI from region II, 
destroys the unstable hmit cycle. 

Region II contains the same fixed points as region VI but no limit cycle. 

The saddle node bifurcation line, which separates region VI from region IV, 
destroys an unstable fixed point and the saddle point. 

The subcritical Hopf bifurcation line, which separates region VI from region 
VII, destroys the unstable limit cycle and turns a stable fixed point unstable. 

Region VII contains only four fixed points (one saddle, two unstable and one 
stable fixed point). 

The saddle node bifurcation line, which separates region VII from region III, 
destroys an unstable fixed point and the saddle point. 

The saddle node bifurcation line, which separates region V from region VIII, 
generates a stable fixed point and a saddle point. 

Region VIII contains four fixed points (one saddle point, one unstable and two 
stable fixed points) and two limit cycles (one stable and one unstable) 
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- The homoclinic bifurcation line, which separates region VIII from region VI, 
destroys the stable limit cycle. 

- The supercritical Hopf bifurcation line, which separates region VIII from re- 
gion IX, destroys the stable limit cycle and turns an unstable fixed point 
stable. 

• Region IX contains four fixed points (one saddle and three stable fixed points) 
and one unstable limit cycle. 

- The saddle node bifurcation line, which separates region IX from region IV, 
destroys a stable fixed point and the saddle point. 

- The homoclinic bifurcation line, which separates region X from region VI, 
generates a stable limit cycle. 

• Region X contains the same fixed points as region VI but two limit cycles (one 
stable and one unstable). 

- The saddle node bifurcation line for limit cycles, which separates region X 
from region II destroys both limit cycles contained in region X. 

Phase portraits for typical parameter values in the regions II and V-IX are displayed 

in Fig. 7. <Fig.7 

Our bifurcation diagram has an intricate structure especially as far as the two 
homoclinic bifurcation lines AiSi and A2S2 are concerned. The insert of Fig. 6 which 
resolves these curves on a finer scale shows that they end up in different points 5*1 
and 5*2. This is a consequence of the fact that the two bifurcations in question play 
essentially different rolls with respect to the global phase portrait of the system 
(cf. Fig. 8). Fig. 8a shows a phase portrait on the line AiSi. One realizes that the < Fig. 8 
homoclinic orbit encloses one unstable fixed point. Fig. 8b shows the same phase 
portrait on the other homoclinic bifurcation line ^25*2. In this case the homoclinic 
orbit encircles two fixed points. Both situations cannot be deformed smoothly into 
each other and are therefore topological distinct^. We add that in Figs. 8a and 8b 
the stable manifold of the saddle point may be thought of as employing just its two 
different branches of the unstable manifold in building the two homoclinic orbits. 

In addition our bifurcation diagram reflects the topology of the phase space. 
For example region IV is bounded by two Hopf bifurcation lines. The upper line 
generates an unstable limit cycle. It is destroyed if one crosses the lower subcritical 
Hopf bifurcation line. As both Hopf lines are connected to different Arnold-Takens- 
Bogdanov points the fixed points involved in the Hopf bifurcations are different. 

^Of course one has to be careful when speaking about the interior of a curve because our phase 
space is a sphere and not the plane. In the present case we define the interior of the homoclinic 
orbit as those part of the phase space which contains the angle between the stable and the unstable 
eigendirection of the saddle point. 
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Hence the limit cycle wanders continuously from one fixed point to the other. Such 
a scenario is impossible if the phase space is a plane, e.g. for a mechanical oscillator. 

If we change the third bifurcation parameter h± the bifurcation diagram changes 
smoothly. However some qualitative modifications occur. Using the same notations 
as in the previous analysis we describe the major differences. Increasing h± the 
Hopf bifurcation line which begins at A2 flattens and finally divides region VIII and 
IX into two parts (cf. Fig. 9). In contrast to regions Villa and IXa regions Vlllb <Fig.9 
and IXb contain an additional unstable limit cycle. At slightly higher values of the 
driving field h± the same phenomenon occurs in region V (cf. Fig. 10). <Fig.lO 

A dramatic change in the bifurcation diagram occurs at the Arnold-Takens- 
Bogdanov point A2 by increasing the driving field further. Fig. 11 shows the diagram <Fig.ll 
for h± = 0.8. The major differences to the previous case (cf. Fig. 10) are 

• The Hopf bifurcation line starting at the point A2 has become supercritical. 

• The homoclinic bifurcation line lies above the Hopf bifurcation line. 

• The saddle node bifurcation for limit cycles now ends up in a degenerated 
Hopf bifurcation point instead of a degenerated homoclinic bifurcation point. 

Hence the Arnold-Takens-Bogdanov bifurcation has changed its type. In terms of 
the normal form 



this bifurcation can be attributed to a change in sign of the coefficient b. This 
degenerated Arnold-Takens-Bogdanov bifurcation is discussed in the mathematical 
literature on an abstract level |T^. In order to keep the discussion self contained as 
well as accessible to physicists we have developed a much more elementary approach 
which nevertheless yields the complete bifurcation diagram of this codimension three 
bifurcation point. It is given in the appendix. 

An additional bifurcation of (at least) codimension three occurs a.t h± = 1. If 
we approach this point the saddle node bifurcation lines collide at 7 = 0, A = and 

generate two new cusp points C3 which initially coincide (cf. Fig. 12). The center <Fig.l2 
manifold of this bifurcation is one dimensional because only a single eigenvalue 
vanishes. On a further increase of the driving field the cusp points C3 separate. 

Fig. 13 shows the bifurcation diagram. <Fig.l3 

4 Conclusion 

Our analysis of a driven spin system has covered both aspects of the dynamics, the 
Hamiltonian as well as the dissipative case. The undamped linearly polarized driven 



X = y 

y = ax^ + hx y 



(49) 
(50) 
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system discussed in section 2 has shown the typical behaviour of a nonintegrable 
Hamiltonian system. Our analysis has revealed that only a finite number of reso- 
nances occurs in each order of the perturbation expansion. This is a consequence 
of the direct coupling of the driving field to the action angle variables (cf. eq. (H)) 
which are in a certain sense the natural variables of the dynamical system. In this 
respect our system differs from e.g. mechanical oscillators where higher harmonics 
come into play even at low order resonances. The local phase space structure near 
these resonances including their width was described by a pendulum equation. In 
addition the phase portrait has been analysed using a global perturbative approach. 
It is in good agreement with numerical simulations. 

The investigation of the dissipative system subjected to a circularly polarized 
driving field has revealed a surprisingly rich bifurcation scenario. Beside local (sad- 
dle node and Hopf ) and global bifurcations (homoclinic and saddle node bifurcations 
of limit cycles) of codimension one we have found five different bifurcations of codi- 
mension two and even two bifurcations of codimension three. From our analysis, 
which has been performed to a large extent analytically, we draw the conclusion that 
our bifurcation diagram is complete. We notice that for a closely related model, i.e. 
an antiferromagnet made up of two homogeneous sublattices, a partial bifurcation 
diagram has been established in pO| . 

It is widely believed and supported by our results that the Hamiltonian and the 
dissipative dynamics are entirely different. But the limit of weak damping and the 
emergence of dissipative dynamics from the Hamiltonian phase space structure is 
poorly understood. The system analysed in this article provides a physical model 
to study this issue. 

Our investigations constitute a necessary prerequisite for a systematic study 
of spatially inhomogeneous states of the full Landau-Lifshitz-Gilbert equation. A 
rather complex behaviour even under the very restrictive assumption of homoge- 
neous magnetization was found. We are convinced that the full spatially inhomo- 
geneous equation, including exchange as well as dipolar interaction, provides the 
correct theoretical description of strongly driven magnetic systems far from equilib- 
rium. Work in this direction is in progress and will be published elsewhere. 
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Appendix: The degenerated Arnold— Takens—Bog- 
danov Bifurcation 



We present here an elementary discussion of the degenerated Arnold-Takens-Bog- 
danov bifurcation. The mathematical background for this codimension three bifur- 



cation was worked out in ref . |19[ . Hence we skip the formal proof of the universality 
of the unfolding but concentrate on the computation of the bifurcation manifolds in 
the full three dimensional parameter space. 

We consider the case that in the well known normal form of the Arnold-Takens- 
Bogdanov bifurcation ( ^ , |50| ) a degeneracy of the second order terms occurs, which 
means that the coefficient b vanishes. In order to establish the normal form of this 
higher order codimension point we start from the general equations of motion 



X2 



/i(x) 

/2(X) 



oo J 



+ EE 



(lijX^X2 



U-i) 



j=3 i=0 
oo j 



X 



•ijX-^X2 



U-i) 



(51) 
(52) 



j=3 i=0 



which coincide up to second order with the degenerated situation just described. 
They contain all possible contributions beyond the second order. In the spirit of 
normal form calculations we are looking for a coordinate transformation x = h(y), 
h(0) = including a rescaling of the time t ^ (t such that the system f^^2[) 



becomes "as simple as possible" |]16j. In the new coordinates we get 

y = g(y) = C(^h(y))'^f(h(y)) . 

Substituting the expansions 



(53) 



hi = cii^/i + E E Cu?/i?/2 * 

j=2 i=0 
4 j 

h2 = duy2 + E E ^iiViyi'' 

j=2 i=0 



(54) 



into equation (p3D and keeping terms up to fourth order we determine the coefficients 
Cij and dij in such a way that as many terms as possible vanish in equation (|53| ). 
The necessary algebraic manipulations are quite extensive. They were performed 
with the help of a computer program applying symbolic mathematics We end 
up with the normal form 



yi 



y2 



o{\: 



yf + a yt + yty2 + 0[\y 



(55) 
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where a is a constant of order unity which remains undetermined. 

In extension and analogy to the two parameter unfoldings of the standard Arnold- 



Takens-Bogdanov bifurcation [Il6| we now introduce the following three parameter 
unfolding 

iji = z/2 + o(|y|') (56) 

2/2 = AH +/i2 2/2 + /U3 Z/1 ?/2 + 2/? + a + 2/1 2/2 + (|y|^) • (57) 



The bifurcation manifolds in the three dimensional parameter space (/xi, /i2, /^s) will 
be computed in the sequel. 

We begin our discussion with an analysis of the local codimension one bifurca- 
tions. Following the lines of section 3.2 we obtain from eqs.(^7|) easily the saddle 
node bifurcation manifold 

/if = . (58) 
Evaluating eqs.(^) a little algebra yields the Hopf bifurcation manifold 

/i2°'''(/ii,/i3) = v/^/^(At3 - /^l) + (/^3 /il,/^?) ■ (59) 
Fig. 14 shows both manifolds in the three dimensional parameter space. < Fig. 14 

To study the global bifurcations a rescaling is performed which differs from 
the blowing up used in the treatment of the ordinary Arnold-Takens-Bogdanov 
bifurcation [|16]. One has to choose 



yi = u, y2 = e^ V, t = e t, (60) 



in order to ensure that the rescaled system becomes Hamiltonian in lowest order 
and that all non Hamiltonian contributions scale with the same order of magnitude 



of e. We obtain from eqs. (|56| , |57|) 



u = V 

V 



+ O (e^) (61) 
vi + u^ + au'^ + V2V + 6^ v^uv + V + [e^^ . (62) 



where the dot now denotes the derivative with respect to the new time r. We note 
that the term proportional to is Hamiltonian. We are now dealing with a four 
parameter problem where the parameters i^i, 1^2, ^^3 are of order 0{1) and e is small. 
In the limit e — eqs. (|6l|j6^) yield an integrable Hamiltonian system with the 
Hamiltonian 



H{u, v) = — -uiu-— . (63) 



Fig. 15 shows the corresponding phase portrait. It exhibits closed orbits and a ho- < Fig. 15 



moclinic loop Fq with the energy H{u,v) = ^J—uf. 
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The singular transformation (pOf ) has blown up the degenerated fixed point into 
a Hamiltonian system because the limit e —>■ implies fii, fi2, f^s ^ 0. By perturbing 
the solutions of the Hamiltonian system (|63D we are able to uncover the behaviour 



of eg s .( 1611 , 1621 ) for parameter values /ii, /i2, /U3 close to the origin. 

Homoclinic bifurcations occur at those parameter values z/i, 1/2, 1^3 for which the 
homoclinic orbit of the Hamiltonian system persists if the perturbation is switched 
on. This problem can be treated analytically by Melnikovs method (e.g. ref . |ll6| ) . 
The (nondegenerated) zeroes of the time independent Melnikov function 

M(z/i, z/2, Z/3) = / (h'2V + PsUV + v) du 

8(/-4z/f 



385 V 309V-z/? + 231z/2-165z/3v^j (64) 

signal the occurrence of the homoclinic bifurcation. In view of the scaling ( |60|) this 
condition reads 



he _ lOSy-Zif + 55 /i3 

The surfaces at which the Hopf, saddle node and homoclinic bifurcations take place 

are depicted in Fig. 16. < Fig. 16 

Finally we show that a further global bifurcation, a saddle node bifurcation 
of limit cycles, comes into existence for the degenerated Arnold-Takens-Bogdanov 
bifurcation. The subharmonic Melnikov function 

Ms(i'i,i'2,^3,E) = / (h'2V + I'suv + v) du (66) 



provides an appropriate tool for studying this issue |16|. Here V{E^vx) denotes 



a trajectory of the Hamiltonian system (|63|) with energy E. If this function has a 
nondegenerated zero for a certain value of the energy then the corresponding periodic 
orbit persists in the perturbed system (I6II J64 ). Moreover a doubly degenerated 
zero of Ms indicates a saddle node bifurcation of limit cycles at the corresponding 
parameter value. Hence the saddle node bifurcation manifold is determined by the 
equations 

= M4Z/1, Z/2, z/3, i?) = V2 / , du (67) 



E + z/im + mV3 



Here the limits of the integrals are given by the first two zeroes of the argument of 
the square root: 

u^ 1 

E + viu + — = -{u — a){b — u){c — u) a<b<c . (69) 
o o 
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We solve the linear system (|67| , |68D for the two parameters z/2 and z/3 and obtain 
a parametric representation of the saddle node bifurcation surface of the form 

(Z/2,Z/3) = F(Z/1,E). 

Fig. 17 contains all bifurcation manifolds in a single diagram. We draw attention < Fig. 17 
to the fact that the saddle node bifurcation manifold for limit cycles meets the Hopf 
surface in a line of degenerated (local) Hopf bifurcations, whereas its contact with 
the homoclinic surface yields a line of degenerated (global) homoclinic bifurcations. 
Both lines represent codimension two bifurcation manifolds. Joyal |Q describes a 



related phenomenon in a different context. In addition we recall that the bifurca- 
tion scenario described at the end of section 3.3 is recovered if one considers two 
dimensional cross sections (e.g. at constant 7x3) of our three dimensional bifurcation 
diagram. In closing we restate that the bifurcation set presented in Fig. 17 yields 
the complete bifurcation diagram of the degenerated Arnold-Takens-Bogdanov bi- 
furcation. 
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Figure captions 



Fig.l Phase portrait of the Hamiltonian at parameter values eh± = 0.01, a = 1, 
= 0.5, and oo = 0.2. 

Fig. 2 Level lines of the approximate second order invariant / at parameter values 
eh± = 0.01, a = 1, Hz = 0.5, cu = 0.2 and at time t = 0. 

Fig. 3 Level lines of the approximate second order invariant / at parameter values 
eh_i = 0.025, a = l, Hz = 0.5, u = 0.2 and at time t = 0. 

Fig. 4 Poincare map at time t = of the full Hamiltonian system @ at parameter 
values eh± = 0.01, a = 1, = 0.5 and uj = 0.2. 

Fig. 5 Poincare map at time t = of the full Hamiltonian system at parameter 
values eh± = 0.01, a = 1, = 0.5 and uj = 0.2. 

Fig. 6 Bifurcation diagram of eq.f|37f) at parameter values h± = 0.1 and F = 0.1. 
The bifurcation lines are denoted as follows: sn: saddle node bifurcation. Ho: 
supercritical Hopf bifurcation, ho: subcritical Hopf bifurcation, he: homo- 
clinic bifurcation, and si: saddle node bifurcation for limit cycles. Codimen- 
sion two points are labeled by capital letters: Ai: Arnold-Takens-Bogdanov 
bifurcation, B: degenerated homochnic bifurcation, Cf. Cusp bifurcation, 
H: degenerated Hopf bifurcation, and Sf. saddle node connection. The 
numbers refer to the regions described in the text. The region (A, 7) G 
[—0.7,-0.69] X [0.39,0.395] is magnified in the small insert. Note that to 
the left hand side of the point B there exists a line of a homoclinic bifurca- 
tions ending up in point 5*2 and a line of a saddle node bifurcation of limit 
cycles ending up in point H. 

Fig. 7 Phase portraits in stereographic projection for several points in the bifur- 
cation diagram Fig.6. The numbers refer to the regions in the bifurcation 
diagram. The phase portraits show actual solutions of the differential equa- 
tions. Parameter values (A, 7) are as follows: (II) (—0.1,0.1) (V) (0.75,0.5), 
(VI) (-0.1,0.6), (VII) (0.5,0.9), (VIII) (-0.65,0.6), and (IX) (-0.65,1). In 
all cases there exists an additional fixed point on the lower hemisphere which 
is not displayed. 

Fig. 8 Phase portraits at the two different homoclinic lines (a) AiSi and (b) A2S2 
(cf. Fig.6). Parameter values are chosen as (a) A = —0.6, 7 = 0.604 and 
(b) A = —0.4, 7 = 0.4009. An additional fixed point exists on the lower 
hemisphere which is not displayed. 

Fig. 9 Bifurcation diagram of eq.(l37|) at parameter values h± = 0.2 and F = 0.1. The 
notation is the same as in Fig.6. 
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Fig. 10 Bifurcation diagram of eq. (|37|) at parameter values /i^ = 0.5 and F = 0.1. The 
notation is the same as in Fig.6. 



Fig. 11 Bifurcation diagram of eg . (|37|) at parameter values h± = 0.8 and F = 0.1. The 
notation is the same as in Fig.6. The lower part shows a magnification of the 
region near the codimension two point A2. 



Fig. 12 Bifurcation diagram of eq.(|37[) at parameter values h± = 1 and F = 0.1. The 
notation is the same as in Fig.6. The insert shows a magnification of the region 
(A, 7) G [-0.02,0.02] X [0.2,0.4] 



Fig. 13 Bifurcation diagram of eq.(|37D at parameter values h± = 1.1 and F = 0.1. 
The notation is the same as in Fig.6. The insert shows a magnification of the 
region (A, 7) G [-0.0015,0.0015] x [0.415,0.445] 



Fig. 14 Partial bifurcation set of system ( 156|J57| ) containing Hopf and saddle node (sn) 
bifurcation manifolds. 



Fig. 15 The phase portrait of the Hamiltonian (163| ) at z^i = — 1. Fq denotes the 
homoclinic loop. 



Fig. 16 Partial bifurcation set of system ( |56|, |57[) containing the Hopf, saddle node (sn) 
and homoclinic (he) bifurcation manifolds (cf. Fig. 14). 



Fig. 17 Complete bifurcation set of system(p6ip7[) containing the Hopf, saddle node 
(sn) and homoclinic (he) bifurcation manifolds as well as the saddle node 
bifurcation manifold of limit cycles (si) (cf. Fig. 16). 
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